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Abstract 

A Monte Carlo code has been developed to simulate the operation of the planned 
detectors in OMNIS, a supernova neutrino observatory. OMNIS will detect neutrinos 
originating from a core collapse supernova by the detection of spalled neutrons from 
Pb- or Fe-nuclei. This might be accomplished using Gd-loaded liquid scintillator. 
Results for the optimum configuration for such modules with respect to both neutron 
detection efficiency and cost efficiency are presented. Careful consideration has been 
given to the expected levels of radioactive backgrounds and their effects. The results 
show that the amount of data to be processed by a software trigger can be reduced 
to the < 10 kHz region and a neutron, once produced in the detector, can be 
detected and identified with an efficiency of > 30%. 

Key words: Supernovae; Monte Carlo; neutrino detector; neutron detector; liquid 
scintillator 

PACS: 14.60.Pq, 29.40.Mc, 29.85.+C, 97.60.Bw 



1 Introduction 



The electron-antineutrinos observed by the Kamiokande (H) and 1MB (g) 
laboratories from supernova 1987-A confirmed that neutrinos play an impor- 
tant role in core collapse supernovae. The vast majority of the gravitational 
binding energy of the protoneutron star is radiated away by neutrinos. The 
observed v e were detected over a few seconds, much longer than the weak 
interaction time-scales which are thought to be required to produce them. 
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This confirmed the main features of the standard model of core-collapse su- 
pernovae (Q) and suggests that the behavior of the neutrino opacity in matter 
beyond nuclear densities is understood (ffl). As has been shown in numerous su- 
pernova simulations (see, e.g. (|)), the shock that will produce the supernova 
stalls, thus disfavoring a prompt mechanism (|6|). Although there are many 
uncertainties in the calculations, indications that the shock has to be revived 
by neutrinos are strong (0). 

An increasing body of data from the atmospheric neutrino experiments 
at Super-Kamiokande (S), 1MB (Rl), Soudan2 (|T0D and Kamiokande (|TT|), ex- 



hibit the disappearance of //-neutrinos ([12]) . Solar neutrino experiments such 
as Homestake (|T|), Kamiokande (0), GALLEX flTJ), SAGE (0) and Super- 
Kamiokande (|I"7D also suggest a flux of u^s inconsistent with the Bahcall- 
Pinsonneault standard solar model. The most favored solution to these dis- 
crepancies is an extension of the standard model of particle physics to one 
in which the neutrinos have mass and are allowed to oscillate between flavors 
via the vacuum and/or the Mikheyev-Smirnov-Wolfenstein (MSW-) oscillation 
mechanisms flTB|). Direct measurements of the neutrino masses (|H]) yield upper 
limits of m Ve < 3 eV (90%CL), m Vlt < 170 keV (90%CL) and m Ur < 18.2 MeV 
(95%CL). However, the tremendous neutrino luminosities of supernovae pro- 
vide a unique scenario for extremely long baseline measurements, allowing 
a vastly increased region of neutrino mixing parameter space to be probed, 
and placing more stringent limits than can otherwise be achieved on neutrino 
masses and mixing angles. 

The astrophysical r-process may also depend critically on the properties 
of neutrinos. The leading candidate for its site is the low-density, high-entropy 
bubble of a core-collapse supernova, created in the trail of the outgoing shock- 
wave by material lifted off the surface of the protoneutron star and heated by 



neutrinos ^Tj). The details of the interaction between the neutrinos and 
matter in the bubble are highly nonlinear, and may require several types of 
MSW-type oscillations (El [23]). 



2 Motivation and Properties of OMNIS 



2. 1 The Role of OMNIS in the Supernova Neutrino Observatory Community 



An observation of the fluxes of different neutrino flavors is vital for the 
understanding of core-collapse supernovae and helpful for the investigation 
of new physics beyond the standard model. In this context, the need for a 
large neutral-current neutrino detector has become increasingly obvious in 
recent years (^|). Accordingly, OMNIS is being designed to provide a much 
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larger sample of /x- and r-neutrino events than would be provided by other 
detectors, allowing a significantly more detailed measurement of the neutrino 
spectra for different flavors, resolved in time, and the stellar conditions that 
produced them. Given the planned quantities of iron and lead, and with the 
detection strategy suggested here, a typical supernova at the center of the 
Galaxy would lead to about 2000 observed events in OMNIS ( fffi) (assuming no 
oscillations), predominantly from neutral- current reactions. These events will 
provide a window for looking well beyond the peak of the supernova neutrino 
burst to study the evolution at times up to pa 10 seconds after the start of 
the pulse, a potentially important diagnostic of the cooling mechanism. A few 
thousand events are also adequate for accurate neutrino mass measurements. 

Most present neutrino detectors, such as Super-Kamiokande, LVD, MACRO 
and AMANDA, are primarily sensitive to charged-current channels, involv- 
ing mainly the reaction of i7 e with protons in water or mineral oil. Super- 
Kamiokande will also see some e-neutrinos and some /i- and r-neutrinos through 
the detection of 7-rays produced in the neutral-current channels 1& 0{y, i/p/wy), 
as suggested by Langanke and Vogel fl2"6p. The total number of detected 
neutral-current events due to \i- and r-neutrinos for a core-collapse supernova 
at a distance of 8 kpc has been estimated to be around 300 or 560, depending 
upon the temperature and chemical potential in the neutrino Fermi-Dirac dis- 
tribution assumed. Recent changes to the data acquisition that result in a lower 
detection threshold might improve this to around 1000 events (f?7|). However, 
this yield comes from a single reaction threshold and so the unknown neutrino 
energy distribution cannot be inferred. SNO, which is designed around a tank 
of heavy water, would be expected to produce around 750 break-up events of 
deuterium following neutrino neutral-current scattering for a supernova at the 
same distance (p8|). It is on the other hand scheduled for a lifetime at the low 
edge of the mean time between Galactic supernovae (10-30 years, (|3CD), even 
though SNO's scheduled time has been extended from two to ten years. 

It has been shown (|28| ) that, based on time-of-flight information, Super- 
Kamiokande and SNO will have resolutions of m v , w 50 eV and m v , « 
30 eV, respectively, although again the improved statistics due to a lower 
detection threshold will improve considerably the value for Super-Kamiokande. 
Due to both the greater number of neutral-current events detected, and to its 
superior timing capability, OMNIS is expected to have a still better mass 
resolution of around 20 eV. 

A particularly interesting scenario is that in which a Galactic core-collapse 
supernova proceeds to a black hole. This should produce an abrupt termina- 
tion of neutrino flux as well as several remarkable measurements. Firstly, the 
signal preceding the cut-off, and when the signal occurs may depend on the 
protoneutron star evolution, and hence the equation of state (p9[). Secondly, 
as the black hole expands outward from the center of the star it will envelope 
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successive neutrinospheres fl29|). One might therefore expect the z/ M / T -flux to 
be terminated earlier than the z7 e -fTux which would in turn occur before the 
termination of the i/ e -flux. The difference in time between these terminations 
has been estimated to be of order 1 ms (|29l). If the supernova is close enough, 
OMNIS would be able to observe these differences, and thus could infer the 
structure of the neutrinospheres. Lastly, the termination of the luminosity 
provides a well defined feature that would provide an opportunity for time- 
of-flight mass measurements of the neutrino mass with unparalleled accuracy, 
down to limits of 2 eV for z/ e in Super-Kamiokande and as low as 4 eV for fi- 
and r-neutrinos from OMNIS (|30|). 



2.2 Expected Neutrino Signal from a Core-Collapse Supernova 



In order to make estimates of the response of the OMNIS detectors to 
a supernova, we have included standard expectations for the incident neu- 
trino intensities, energy distributions and interaction cross sections. The bulk 
of the neutrino luminosity from a core-collapse supernova is emitted during 
the Kelvin-Helmholtz cooling phase, which is marked by an exponentially de- 
caying neutrino flux with a time constant of ~3 seconds fl3T|). Each flavor is 
described by a Fermi-Dirac distribution with a chemical potential at or near 
zero. At the relevant energies, the /i- and r-neutrinos only interact via neutral- 
current processes, and therefore decouple deeper within the core (at a smaller 
neutrinosphere radius, the term being defined as the radius from the center at 
which neutrinos decouple from the dense matter), leaving them with the high- 
est energies, estimated to be around < E v . >m 25 MeV or T « 8 MeV. 
Electron-neutrinos decouple at a larger neutrinosphere radius, with the net 
numerical superiority of neutrons to protons resulting in the e-antineutrinos 
having higher energy than the e-neutrinos (< E Vg >~ 16 MeV or T « 5 MeV 
versus < E Ue >^ 11 MeV or T « 3.5 MeV (||)). 



2.3 Neutrino Cross Sections Relevant to OMNIS 



OMNIS will consist of several ^ kT modules of lead and iron, which will 
serve as target materials. Some of the incoming neutrinos from a supernova 
will undergo charged- current or neutral-current interactions with the compos- 
ite nuclei, resulting in neutron emission. The cross sections for such reactions 
depend on the neutrino energy, the Q-value for the reaction and the avail- 
ability of suitable energy levels in the target nucleus. There are no available 
experimental data for the cross sections governing these reactions, with the ex- 
ception of that for 56 Fe(V e , e~) 56 Co measured by the KARMEN collaboration 



but theoretical studies have been conducted by two groups ([34]; |35|). De- 
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Table 1 

Comparison for single- and double-neutron events liberated from iron and lead for 
the different reaction channels (per kT) 



Material, Event Type 


CC-i/ e 


CC-Pe 


NC-u e 


NC-P e 


NC-i/j 


total 


Pb, single-n, no osc. 


59 





8 


37 


677 


781 


Pb, double-n, no osc. 


26 








1 


20 


47 


Fe, single-n, no osc. 


4 


5 


2 


6 


146 


163 



spite qualitative agreement between the two groups, some uncertainty remains, 
due to the largely unknown location of the Gamow- Teller resonance relative 
to the neutron separation energies. These separation energies are 10.67 MeV 
for 56 Fe and 7.19 MeV for 208 Pb for single-neutron emission and 19.78 MeV 
and 13.31 MeV, respectively, for double-neutron emission. The single-neutron 
channels in lead and in iron and the two-neutron channel in lead constitute 
three observables closely related to the neutral-current thresholds that can be 
used to deduce the relevant contributions of the various neutrino flavors via 
the different energy dependencies of their cross sections. Both SNO and Super- 
Kamiokande will provide additional reaction thresholds and their data could 
be used in collaboration with those from OMNIS to map out the neutrino 
distributions. 



2.4 Expected Number of Events in OMNIS 

All event numbers presented in this work are based on a "standard" su- 
pernova in the center of the Galaxy (distance = 8 kpc, E to t = 3 x 10 53 ergs). 
The number of neutrons liberated per kT in OMNIS modules of lead or iron 
is indicated in table 1. The cross sections for lead were taken from fl3"4"D and 
for iron from fl35|) . 

2. 5 Influence of Neutrino Oscillations on the OMNIS Yields 

One possible scenario is the oscillation between \i- or r-neutrinos and 
electron-neutrinos. The resulting higher energy for the latter, combined with 
the large charged- current cross section, leads to a significant increase in the 
overall yield in both the iron and the lead and to a dramatically increased 
number of double-neutron events in lead (see also table 2, which shows the 
result for an oscillation of all /x-neutrinos into e-neutrinos and all e-neutrinos 
into /i-neutrinos) . The latter channel therefore serves as a sensitive neutrino 
thermometer and/or indicator for certain types of neutrino oscillations. Fig- 
ure 1 shows the expected number of liberated neutrons per kT as a function of 
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Table 2 

Comparison for single- and double-neutron events liberated from iron and lead for 
the different reaction channels (per kT) with full oscillation i/„ u e . 



Material, Event Type 


CC-v e 


CC-v e 


NC-z/ e 


NC-P e 


NC-u x 


total 


Pb, single-n 


826 





184 


35 


516 


1563 


Pb, double-n 


1852 





6 


1 


15 


1874 


Fe, single-n 


57 


5 


37 


6 


112 


217 
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Fig. 1. Neutron Yield per kT Versus Oscillation Probability i>„ <-> v e . 

P, the chance that a neutrino emitted from the supernova as a /i-neutrino has 
oscillated into an electron-neutrino and vice versa. The ability of OMNIS to 
discriminate between possible oscillation scenarios is illustrated in fig. 2, where 
the difference between the number of single-neutron emission events and the 
number of two-neutron emission events observed in lead, is plotted against the 
number of neutron emission events in iron. This demonstrates the additional 
sensitivity that OMNIS has to oscillations involving sterile neutrinos. 



3 Features of the OMNIS detector 



3.1 The OMNISita Test Facility 



OMNIS will consist of approximately 4 kT of iron and 8 kT of lead 
modules. It will be built in the Waste Isolation Pilot Plant (WIPP) near 
Carlsbad, NM, which has been shown, given its background radiation and 
available infrastructure, to be a highly desirable site Q3"6D . At the time of this 
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Fig. 2. Single-Neutron Yields in the Pb and Fe modules (each per kT) for different 
Oscillation Schemes. 

paper, a test module, OMNISita fl37D, is being constructed at the proposed 
site in order to test the predictions of the Monte Carlo code presented here, to 
have a test facility for the components of the OMNIS detector, and for further 
measurements of the radioactive and cosmic ray backgrounds. 

3.2 Gd loaded Liquid Scintillator for Neutron Detection 

Once neutrons are emitted, they interact with the neutron source material 
(Fe or Pb) and the neutron detector, assumed in the present study to be 
organic liquid scintillator. This serves three functions: neutron moderator, sink 
and detector. The use of Gd in liquid organic scintillators is well established 
(|38|), it has the highest known neutron capture cross section (|39|) and its 
solubility is high enough to make it the preferred neutron sink at this time. 
The possible alternatives Li and B have not been found to be soluble in organic 
liquids in sufficiently high concentrations and there are no known organic 
molecules containing Li or B that have a light yield similar to Gd loaded 
scintillator. 

3.3 Particle Propagation and Event Identification 

The neutrons produced by the neutrino interactions are expected to have 
initial energies of ~ 1 MeV fl35D, hence moderation will primarily occur 
through elastic scattering off the protons in the scintillator. Our calculations 
show that the neutrons lose nearly all of their kinetic energy within oc 100 ns 
(see fig. 3) after emission, resulting in the production of a significant number 
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moderation time (nanoseconds) 

Fig. 3. Neutron Moderation Time-scale 

of scintillation photons. With a 0.1% loading of Gd in the scintillator, capture 
occurs after « 30 /is, see also figure 4, predominantly on 155 ^ 157 Gd. Figures 3 
and 4 were generated with the optimized design for the lead modules presented 
below, but these time-scales are not very dependent on the detector material 
and geometric configuration used. Further, the capture time scales are con- 
sistent with a recent study performed using a prototype scintillator 
vessel and with similar Monte Carlo studies conducted in the past (|41| ; |42|) . To 
enable discrimination against cosmic rays and 7-rays produced by the decay 
of radioactive impurities in the lead or iron or other detector components, 
a double-pulse technique with a time window between prompt and capture 
signals of ~ 50 fis will be used. The characteristic second pulse, caused by 
the emission of 7-rays following capture on Gd or H in the scintillator vessels, 
is used to identify the neutron pulse. The time of neutron liberation is then 
known within a few 100 ns. Although pulse shape discrimination might be 
considered as a means to identify the neutron induced events, this would be 
impractical for detectors as large as those required for a facility of the size of 
OMNIS. 



3.4 Design Constraints 

Generally, the more fine-grained neutron sources and detection elements 
are, the higher will be the neutron detection efficiency, since neutron energy 
losses and captures before they reach the detector will be minimized (|3"8D . 
One important cost overhead though will be the data acquisition system; this 
imposes a practical upper limit of several hundred independent scintillator 
units per kT. It is important for the lead to be free of antimony, which acts as 
a significant neutron sink in concentrations above ~ 1%. Structural stability 
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Fig. 4. Neutron Capture Time-scale 

dictates lead walls with thicknesses of at least ~ 10 cm, interlaced by layers of 
scintillator with a thickness determined by the neutron attenuation length in 
organic liquids, which is several cm. Rectangular acryllic tubes are assumed 
as scintillator containment vessels, which guide the scintillation light via total 
internal reflection off the acryllic/air interface. Their length will be limited by 
the attenuation length for photons in the scintillator, which is about 450 cm 
(|43|). The rectangular shape maximally fills out the scintillator walls, and 
can be manufactured for a considerably lower price than comparable units 
with different geometries and materials. For the regions on the tank ends not 
occupied by the photomultiplier tubes, a highly reflective coating of aluminum 
or T12O is planned. A similar design was used for the Palo Verde (fJ3j) reactor 
experiment, which utilized acryllic vessels with 0.1% Gd loaded pseudocumene 
(PC) as a liquid scintillator in a solution with 60% mineral oil. The present 
analysis has assumed a similar material which has a high proton concentration 
(|j ~ 1-64), making it a good moderator. However, to increase the range over 
which the thickness of the detector walls can vary, stacking of the scintillator 
vessels in single-, double- or triple-width columns has been considered. To 
provide shielding from external radiation sources and reflect escaping neutrons 
back into the detector, each module is encased in an external lead or iron hull. 
Figures 5 and 6 show a schematic design for a \ kT lead module. 



4 DAMOCLES - Model for OMNIS 

4-1 Basics of Monte Carlo Code Development 

The Monte Carlo code for OMNIS was designed to provide: 
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Fig. 5. Front View of an OMNIS Lead Module 

(1) Modularity: The detector consists of three principal regions, the scintilla- 
tor vessels, the lead walls and the outer hull which have to be simulated 
in a large number of different configurations, necessitating easy control of 
the detector configuration. In particular, the code handles the operation 
of a large number of identical cells. 

(2) Parameter Control: The development of OMNIS requires the control of 
all parameters regarding geometry, data acquisition and the materials 
used, including the composition of the scintillator and possible sources 
of radioactive background in all detector components. That includes the 
cost of a given detector configuration, which is necessary to optimize cost 
efficiency 

(3) Data Flow Control: In order to design the best trigger mechanism, ex- 
traction of any kind of information about particle positions and energy 
deposited in the detector was accommodated. 

(4) Portability: Since the simulations were performed on both single- and 
parallel- processor architectures, the code had to be portable among dif- 
ferent platforms. 
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SCINTILLATOR VESSEL 
Fig. 6. Top View of Part of an OMNIS Lead Module 

4-2 Chronology of Events in the Monte Carlo Simulation 



The treatment of events within the Monte Carlo code follows their chronol- 
ogy. Neutrons are first generated with appropriate energies from the spallation 
following neutrino interactions. Each neutron is then tracked through the de- 
tector, taking into account scattering reactions and eventual capture in one of 
the components or its escape from the detector. In the former case, capture 
7-rays are traced through the detector until their eventual capture or escape. 
Neutrons leaving the detector are not tracked any further. The effects of in- 
cluding an additional external reflecting hull were tested beyond that required 
for shielding out external background radiation from the surrounding rock, 
but the resulting changes were negligible. Energy deposited into the scintil- 
lator via neutron or 7-interactions is converted into a number of scintillation 
photons, each of which is then traced until its absorption in the liquid, failure 
to reflect from a surface, or transmission to a photomultiplier tube. 
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4-3 Energy Spectrum of Double-Neutron Events 



In the case of a double-neutron emission event, the available energy is 
distributed to both neutrons Q55D, and phase space arguments tend to equalize 
the energy of both neutrons. A Gaussian energy distribution was therefore 
assumed with the energy for each neutron limited to the total available energy, 
centered at and with a width of half the total energy. The influence of the width 
of the Gaussian on the detection efficiency for double-neutron events in lead 
was investigated. The narrower the Gaussian, the higher the efficiency, since 
the probability for one neutron having very low energy is then small. However, 
the double-neutron efficiency was found not to depend very sensitively on the 
width; the efficiency changed by only 1%, if the width is changed to either 
10% or 200% of the original value. 



4-4 Neutron Interactions 



4-4-1 Neutron Scattering 

The total elastic scattering and absorption cross sections have been ap- 
proximated using polynomials which were fitted to data from the ENDF/B file 



(44). Resonances were averaged over and not treated in detail, since they are 
limited to narrow energy ranges. At the neutron energies with which OMNIS 
will be dealing, ~ 1 MeV, inelastic scattering can be neglected. The expansion 
of the elastic differential scattering cross section in Legendre polynomials can 
be written as 

2/ + 1 

f(cos(8), E) = Z t (—I-a e (E)P t (cos(0))), (1) 



where ai are the (energy dependent) Legendre coefficients. Step functions were 
used based on the known Legendre coefficients as functions of the energy, and 
the series was terminated at £ = 5 for Fe and I = 10 for Pb, beyond which 
contributions are negligible (|4~4"|) . Both C and H (scintillator) exhibit an s-wave 
distribution in the angular dependence of elastic scattering in the center of 
mass system for the energy ranges considered , thus terminating the series 
in eq.(l) after i = 0. Since its concentration is only 0.1%, and its scattering 
cross section is comparable to the host material, we neglect scattering from 
Gd. 
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4-4-2 Neutron Capture 

Neutron capture on Gd is well described by four photons with a total 
energy of 7.937 MeV ( 157 Gd) or 8.536 MeV ( 155 Gd) ([41]) with individual 
energies of between 50 keV and 3.5 MeV (|45|) . Capture on a proton leads 
to the emission of a 2.2 MeV photon. Capture on carbon was negligible, as 
was confirmed by tests. A random 7-ray of up to 4 MeV was assumed to 
be produced by neutron capture on iron or lead. However, our results did not 
change at all if this process was left out completely, due to the high absorption 
cross section in Fe or Pb for 7-rays. For the same reason, any secondary 7-rays 
produced in inelastic scattering off Pb or Fe were found to be negligible. 

4-5 ry-ray Interactions 

The main energy loss channel for 7-rays in the energy range of interest 
here, up to a few MeV, in both the scintillator and the lead or iron, is Compton 
scattering, although all possible effects were included. The cross sections for 
the photoelectric effect, pair production, coherent scattering and Compton 
scattering were obtained from the XCOM database ( |46| ) and fitted for energy 
bins of 10 keV (1 keV for the first 100 keV in the scintillator). For the angular 
dependence in the differential cross section for Compton scattering the Klein- 
Nishina equation for photons scattering off electrons was used (^) . Immediate 
photoelectric absorption was assumed for 7-rays below an energy of 5 keV in 
Fe and Pb and 0.5 keV in the scintillator. Below these respective energies, all 
interaction channels except photoelectric absorption are negligible, and the 
attenuation length is < 0.6 mm in the scintillator, < 0.12 mm in lead and 
< 0.9 mm in iron (flop. 

4-6 Scintillation Photons 

The liquid scintillator used in the Palo Verde experiment had a light yield 
of 56% of that of anthracene ([43"D, the latter being 16 photons per 1 keV (f48"l). 
The neutron energy conversion efficiency due to the higher charge density 
caused by recoil protons is about 40% of that of the 7-ray conversion effi- 
ciency (|49|) . The phototubes were assumed to detect a given photon with an 
efficiency of 20%. Since the tracking of the scintillation photons constitutes 
the bulk of computing time, only 20% of the actual number of photons were 
simulated, but a detection efficiency of 100% was adopted in the phototube. 
Scintillation photons are created with random polarization along either of two 
possible axes. The absorption cross section in the fluid is determined by an 
exponential law with an attenuation length of 450 cm. If the incident angle 
on a surface is larger than the limit for total reflection, the photon is reflected 
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with 100% efficiency. The validity of this assumption was confirmed by intro- 
ducing an additional finite absorption probability when a scintillation photon 
traverses the acryllic wall. Even for a transmission of only 95% for every reflec- 
tion, the efficiency decreases by only 1%. For all other incident angles, every 
reflection off the walls was tested for absorption, with the reflection probabil- 
ity being a function of polarization and incident angle Q47|) . If the photon hits 
one of the scintillator vessel ends, it is either transmitted through the photo- 
multiplier glass face with 98% efficiency or reflected off the highly reflective 
rim surrounding the phototube with an efficiency of 88% for aluminum. The 
latter value was confirmed by measurements conducted at The Ohio State 
University. 



4-7 Event Identification 



In order to discriminate against noise in the photomultiplier tubes and 
the electronics, a threshold was set requiring that one vessel has to register 
at least five photo-electrons on each end within a timebin. The length of a 
timebin was 100 ns, in accordance with the expected resolution of the data 
acquisition system. The threshold of five photoelectrons was chosen because it 
is the lowest pulse that can reliably be detected in a photomultiplier tube, but 
the effect of changing that value was studied as well. If two such pulses occur 
in the same vessel within the time-to-amplitude converter time window, i.e. if 
both ends of one vessel therefore fire twice, an event is registered in that vessel. 
Double-neutron events are identified based on two such events registering in 
separate vessels with the first pulse in each event happening within 100 ns. 



4-8 Intrinsic Background Radiation 



Assuming that only virgin lead will be used in OMNIS, all major internal 
background sources in the lead modules will be elements in the natural decay 
series or their daughters. In lead from the DOE RUN company, a survey of 
samples exhibited background ranging from < 0.001 to ~ 10 ©, with 



most samples towards the lower end of the range. The dominant background 
source is 210 Pb (Ti = 22.2 yr) ([51]). One of its daughters, the short-lived 
210 Bi, emits 1.17 MeV 7-rays Q52T). The absence of significant contributions 
from other isotopes, in particular the 226 Th series, was recently confirmed by 
CEMRC in Carlsbad flSBp. Bremsstrahlung from the various beta-decays can 



be neglected, because its intensity peaks at energies below ~ 0.25 MeV, where 
lead absorbs the photons very efficiently. 

The total radioactivity in iron ranges between 0.00 and 0.17 decays per 
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m Co (1.17 MeV and 1.33 MeV, respectively; 60% of the background) and 40 if 
(1.46 MeV; 30% of the background). An admixture of 10% 226 Th (2.6 MeV 
7-ray from 208 Tl) was assumed for higher energy 7-ray emissions. 

Another background source is the decay of m K in the glass faces of the 
photomultiplier tubes. The available purity ranges from 1 to 100 decays per 
second per kg per face (|55|). 

The detected frequency of false fp events should scale with the back- 
ground rate b as 



where k\ and k 2 are constants depending on the trigger time window and fn 
is the frequency of real neutron events. The linear term in b enhances the 
detection efficiency by providing a missing second pulse to go with either the 
prompt (neutron-) or delayed (7-) pulse. The quadratic term in b describes 
coincident double-pulse triggers from the background, in which both pulses 
within the trigger window are caused by background radiation. It is important 
for the latter term to be small at all times, even when there is no supernova, in 
order to avoid false alarms. During the actual supernova pulse, the background 
contribution is then dominated by the linear term. 

The number of background 7-rays occurring in a time bin, as implemented 
in the code, is Poisson distributed: 



where k denotes the resulting number of background events and A is the num- 
ber of events per time bin. 



5 Efficiency Optimization 

5. 1 Optimization of Detector Dimensions 

The basic philosophy adopted was to optimize the cost efficiency of the 
OMNIS modules rather than the absolute efficiency Cost efficiency here is 
defined as efficiency divided by cost times US% 10 7 , which ranges between 
and 2. The dependence of the detection efficiency and cost efficiency on 



f F = hbf R + k 2 b 2 , 



(2) 



Pk = 




(3) 
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Fig. 7. Efficiency and cost efficiency vs. lead wall thickness for single-neutron events. 

the parameters scintillator vessel dimensions, lead/iron wall thickness, hull 
thickness, photon threshold and neutron energy was investigated. 



5.1.1 The Lead Modules 

The most cost efficient way to detect single-neutron events was found to be 
a design with single scintillator columns with photomultiplier tubes (PMTs) 
of 10 in (25.4 cm) diameter. Double and triple columns of photomultiplier 
tubes with smaller PMTs were also considered. However, scintillator vessels 
with a diameter that is considerably larger than the attenuation length of both 
neutrons and 7-rays dramatically increase the detection efficiency and decrease 
the number of scintillator vessels needed for a given mass of lead. This leads to 
a 50% better cost efficiency for single columns with 25.4 cm-PMTs compared 
to double columns of 12.7 cm-PMTs. For double-neutron events, the greater 
stopping power for larger vessels also leads to a higher cost efficiency, more 
than compensating for the detrimental effect of a less fine-grained detector on 
the detection of two neutrons in two separate detectors. 

The cost efficiency for both single- and double-neutron event detection 
rises as the lead wall thickness increases, then flattens out for a thickness of 
50 cm (fig. 7 and 8), although from a pure efficiency-point of view, somewhat 
thinner lead walls would be preferable. 

As expected from the bulk photon attenuation length in the scintillator, 
a vessel length of 300 cm maximizes both efficiency and cost efficiency for 
single- and especially double-neutron event detection (fig. 9, 10). 

The cost efficiency for single-neutron events might be optimized for ves- 
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Fig. 8. Efficiency and cost efficiency vs. lead wall thickness for double-neutron events. 
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Fig. 9. Efficiency and cost efficiency vs. scintillator vessel length in the OMNIS lead 
module for single-neutron events. 

sels with widths slightly larger than required by the size of the PMTs, the 
fewer number of PMTs necessary might more than compensate for the lower 
efficiency for each single scintillator vessel. There is a slight decrease in the 
cost efficiency with increasing vessel width (by 5 — 10%) and a slight increase 
(by ~ 3%) with vessel height, when increased from 30 cm to 45 cm. However, 
the cost efficiency for two-neutron events requires that the vessel dimensions 
should only exceed the PMT dimensions by a few cm at most, since a decrease 
of the cost efficiency for two-neutron events by ~ 20% is observed when both 
vessel width and height are increased from 30 cm to 45 cm. 

The calculation does not consider the cost for the electronics, although it 
will be small compared to the cost for the scintillator vessels and the PMTs. 
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Fig. 10. Efficiency and cost efficiency vs. scintillator vessel length in the OMNIS 
lead module for double-neutron events. 

However, it decreases with an increasing number of scintillator vessels. We 
have considered using four PMTs with diameters of 12.7 cm instead of one 
25.4 cm-PMT per scintillator vessel. This would provide the same fractional 
coverage of the vessel face area, but on top of the approximately four-fold cost 
for the electronics, one larger PMT is less than four times as expensive as four 
smaller PMTs. 



A thicker lead hull slightly decreases both the detection efficiency and the 
cost efficiency of the detector, since more of the neutrons that are produced 
are lost. We therefore plan to make the hull only as thick as necessary to 
attenuate the external background sources in the WIPP cave. 

The detection efficiency for single-neutron events is constant for neutron 
energies beyond a few hundred keV (fig. 11). For double-neutron events, the 
available energy gets distributed to two neutrons. Hence, the detection effi- 
ciency for those events starts showing an asymptotic behavior only for energies 
above ~ 1 MeV (fig. 12). 

The detection efficiency is still satisfactory, if the pulse threshold has to 
be set higher than the five photoelectrons assumed (fig. 13). However, as ex- 
pected, the double-neutron efficiency decreases more rapidly for higher thresh- 
olds (fig. 14). 



5.1.2 The Iron Modules 

The iron modules have to be optimized only with respect to the cost ef- 
ficiency for single-neutron events. The optimum dimensions of the scintillator 
containment vessels and the wall thicknesses are similar to those for the lead 
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Fig. 11. Efficiency vs. neutron energy in the OMNIS lead module for single-neutron 
events. 
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Fig. 12. Efficiency vs. neutron energy in the OMNIS lead module for double-neutron 
events. 

modules. As in the case of lead, an iron detector with fewer and larger scin- 
tillator vessels is both more efficient and cost efficient. The dependence of the 
cost efficiency on the vessel width and height suggests more elongated vessels 
along the vertical direction than in the case of lead. The detection efficiency 
becomes flat for higher neutron energies than in the case of lead (fig. 15), which 
is due to the fact that more energy is lost in an average elastic collision with 
an iron nucleus. 
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Fig. 14. Efficiency vs. photoelectron threshold in the OMNIS lead module for dou- 
ble-neutron events. 



5.1.3 Summary for Optimized Dimensions 

Table 3 summarizes the optimized dimensions for |A;T lead and iron mod- 
ules, stacked in single scintillator vessel columns (with photomultiplier tubes 
with diameters of 25.4 cm on each end). 



5.2 Influence of the Background 
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Fig. 15. Efficiency vs. neutron energy in the OMNIS iron moduie. 



Table 3 

Summary of optimized detector dimensions for the OMNIS lead and iron modules. 



Material 


wall (cm) 


hull(cm) 


vessel(L(cm)xW(cm)xH(cm)) 


Pb 


50 


15 


300 x 30 x 30 


Fe 


50 


25 


300 x 30 x 45 



5.2.1 The Lead Modules 

The detected false event rate in OMNIS in the absence of emitted neutrons 
is shown in fig. 16. This confirms radiological analyses fl5l|) which show that 
if the bulk impurity rate in lead can be maintained below 0.1 then the 
false event rate in OMNIS will be less than ~ 1 per second. The ability to 
detect neutrino events at late times after the core bounce is determined by 
the influence of the radioactive background on a sparse neutron signal, with a 
typical time interval between two neutrons being long compared to the dead 
time of the detector, taken to be 80 [is. The apparent efficiency of a sparse 
neutron signal relative to the efficiency without background is shown in fig. 
17, which shows that a tolerance of 0.1 is also acceptable in this respect. 

Figures 18 and 19 show the false event rate due to 40 K decays in the pho- 
tomultiplier glass faces without neutrons and the relative apparent efficiency 
in the presence of a sparse neutron signal. A tolerance limit of 100 Bq has 
been adopted for each PMT. 
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Fig. 16. False Event Rate per kT due to Bulk Background in OMNIS Pb Modules. 
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Fig. 17. Apparent Enhancement of Sparse Neutron Signal due to Bulk Background 
in OMNIS Pb Modules. The two horizontal lines indicate a ±5% change relative to 
zero background. 

5.2.2 The Iron Modules 

Due to the higher energy of the background 7-rays and their longer atten- 
uation length in iron, the requirements for its purity have to be more stringent 
than for lead. The resulting tolerance limit for the bulk background in the iron 
used is ~ 0.05 |p (see figures 20 and 21). This also appears to be consistent 
with radiological analyses (|5*^). The allowed background rate in the photo- 
multiplier glass faces is, however, the same as in lead, 100 Bq per PMT (see 
figures 22 and 23). 
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Fig. 18. False Event Rate per kT due to 40 K decays in PMT faces in OMNIS Pb 
Modules. 
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Fig. 19. Apparent Enhancement of Sparse Neutron Signal due to 40 K decays in PMT 
in OMNIS Pb Modules. The two horizontal lines indicate a ±5% change relative to 
zero background. 



5.3 Dependence on Supernova Distance 



If one event strobe is used for an entire detector module, there is a min- 
imum distance of the supernova for which its dead time becomes critical for 
reliable event detection. Assuming one detected neutron event to trigger a 
strobe with 80 /is for a lead module of \ kT (with 60 scintillator vessels), 
dead time losses become significant for supernovae closer than 2 kpc, which 
represents only a few percent of the Galaxy. Table 4, generated with a sim- 
ulated supernova neutrino burst (|56"D, shows the number of detected events 
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Fig. 20. False Event Rate per kT due to Bulk Background in OMNIS Fe Modules. 
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Fig. 21. Apparent Enhancement of Sparse Neutron Signal due to Bulk Background 
in OMNIS Fe Modules. The two horizontal lines indicate a ±5% change relative to 
zero background. 



versus supernova distance. The case without dead time is shown for reference. 
We intend to utilize a sufficiently sophisticated data acquisition system such 
that it can handle the higher neutron liberation rates associated with close 
supernovae in at least one module, thus precluding the excessive loss of high 
resolution data obtainable from a close supernova. 
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Fig. 22. False Event Rate per kT due to 40 K decays in PMT faces in OMNIS Fe 
Modules. 
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Fig. 23. Apparent Enhancement of Sparse Neutron Signal due to 40 K decays in 
PMT faces in OMNIS Fe Modules. The two horizontal lines indicate a ±5% change 
relative to zero background. 



Table 4 



Number of detected neutron events versus supernova distance for 16 -^kT Pb mod- 
ules. 



Dead Time (fis) 


0.20 kpc 


0.50 kpc 


1.0 kpc 


2.0 kpc 


4.0 kpc 


8.0 kpc 


16 kpc 


80 


0.888 x 10 6 


0.293 x 10 6 


95400 


26200 


6740 


1730 


440 





2.86 x 10 6 


0.458 x 10 6 


112000 


27500 


6860 


1740 


440 
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6 Summary 



Once built, OMNIS will fill in the long perceived need for a large neutral- 
current supernova neutrino detector. The neutrino fluxes for the different 
flavors will be measured using the yields from multiple threshold processes, 
single- and double-neutron emission in lead and single-neutron emission in 
iron. The time resolution is only limited by the event rate (and therefore the 
distance) of the supernova. This will reveal unprecedented insight into the 
core-collapse supernova mechanism and might, for the first time ever, open up 
a window onto the black hole formation mechanism, and possibly enable the 
observation of the different neutrinosphere radii for different neutrino flavors. 
In addition, OMNIS is capable of placing new limits on the allowed neutrino 
mass - mixing angle phase space using time-of-flight methods for different 
flavors and the relative yields for the different neutron liberation channels. 

The generic designs for the OMNIS modules are walls of lead or iron, 
interlaced with stacks of vessels containing organic liquid scintillator loaded 
with Gd. A double-pulse trigger with the capture 7-ray pulse giving the con- 
firmation to the preceding neutron pulse in the same scintillator vessel was 
identified as an effective means to discriminate neutron events against the ra- 
dioactive background. Using the Monte Carlo code DAMOCLES, developed 
for OMNIS, optimum configurations for the detector geometry of the iron 
and the lead modules were found with respect to the maximum number of 
events for a given unit cost. The resulting detection efficiencies for | kT mod- 
ules are ~ 38% for single- and ~ 13% for double-neutrons in lead and ~ 34% 
for iron. However, the detection efficiency might improve due to lower possible 
thresholds when off-line analysis of the data is performed. Given the predicted 
numbers of neutrons liberated, we estimate a total of 8 kT of lead and 4 kT 
of iron is needed for the OMNIS detector. 



The false trigger rate caused by the radioactive background in the pres- 
ence or absence of a supernova neutrino pulse was investigated with the result 
that tolerance levels lie within the range found in commercially available ma- 
terials. The criterion was a false rate of ~ — , low enough to still detect a 

sec. ' 

sparse neutron signal. 



The overall neutron trigger rate is below 10 kHz for a | kT module and 
a "standard" Galactic supernova (|56| ; |5T|) , which would enable the use of a 
common trigger gate for the entire detector. If, however, a supernova occurs 
at distances less than ~ 1 kpc, the typical time between two-neutron events 
at the peak of the signal becomes lower than typical dead times of ~ 80 /is, 
in which case common trigger gates for smaller numbers of or even individual 
photomultiplier tubes become desirable in order to measure the peak neutron 
liberation rate, at least for part of OMNIS. 
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